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FOREWORD 


This technical report presents the results of the Laser-Heated Rocket 
Thruster Program. The analysis, design and fabrication were conducted 
by Rocketdyne Division, Rockwell International, during the period 
27 June 1975 to 2 July 1976 as part of National Aeronautics and Space 
Administration, Lewis Research Center, Contract NASS- 19728. The lOkw 
experimental hardware was shipped to NASA-LeRC. 

The NASA-LeRC Project Manager was Mr. S. Cohen and the Rocketdyne Program 
Manager was Dr. S. V. Gunn. 
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SUMMARY 

Conceptually the laser-heated rocket may employ one of two methods of opera- 
tion which are illustrated in Fig. 1 and 2. The first concept has a stationary 
plasma with a recirculating chamber flow. The second concept utilizes a focused 
beam with a flowing plasma and a continuous chamber flow without recirculation. 

The final concept selected in this program was a combination of the two con- 
cepts. The selected concept used a focused beam with a flowing plasma and re- 
circulating flow pattern upstream of the plasma and continuous flow in the 
region of the plasma. 

In this ij-ogram, analysis and design of the 10-kw cxierimental thruster was 
performed in Task I. Physical Sciences, Inc. defined the plasma core size 
and through ^low (stable plasma) and heat transfer influences the 3.45 
X 10 N/M (50 psia) chamber pressure 10-kw thruster was sized for sea 
level operation. The throat diameter of this chamber was 0.1168 cm (0.0468 
inch) and the main hydrogen injector was 8.89 cm (3.5 inches) upstream of 
the throat with a 11.43 cm (4.5 inch) distance from the 4.445 cm (1.75 inch) 
diameter zinc selenide window to the throat. The igniter (plasma arc initia- 
ting electrode mecha'ism) with retracting tungsten electrodes is located 
2.54 cm (1-inch) upstream of the throat. With the predicted peak heat flux 
of 2519 watts/ cm^ (15.4 Btu/in.^sec), the water-cooled chamber would allow 
steady state operation with a 1.7 to 2.3 margin on burnout heat flux. This 
could be accomplished without carbon seed and therefore to avoid the potential 
throat plugging problem, pure hydrogen was selected as the propellant. 

Although the plasma size was predicted with the most sophisticated analytical 
techniques available, a possiblity exists that the plasma size may be sig- 
nificantly larger than predicted. If this occurred in the water-cooled 
chamber, the walls of the chamber would overheat. Therefore, an uncooled 
chamber which could handle a plasma up to four times the predicted diameter f 

was designed. The thrust stand was designed to accommodate cither the water- I 

I 
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Figure 1 . Stationary Plasma With Recirculating 
Chamber Flow Pattern 



Figure 2. Flowing Plasma With Continuous Chamber 
Flow Pattern 



cooled or uncooled chamber using a balance beam concept with a load cell to 
measure thrust. 

Using methods similar to the JANNAF (Joint Array, Navy, NASA, Air Force) 

rocket thrust chamber performance methodology, the predicted specific impulse 

for the 10-kw thruster was 1232. S lb. sec/lb . 

f m 

The 10-kw thruster component designs prepared in Task T were fabricated which 
included the water-cooled chamber, uncooled chamber, injector, igniters, 
and thrust stand. 

Detail design evaluation of several optical train configurations for the 10-kw 
experimental thruster were considered and a pointer cassegrain was selected. 
Beam propagation and chamber window thermal analyses were performed. De- 
tailed drawings of the recommended optical train design were prepared and sent 
to the NASA-LeRC Project Manager. 

The conceptual evaluation of the 5000-kw chamber for space application re- 
vealed that carbon seeding was necessary to reduce the plasma radiation heat 
input. To accomplish this, the radiation absorption layer must have a sig- 
nificantly lower velocity than the primary hot gas. Perhaps the best method 
of introducing the carbon seeded hydrogen gas is to use multi-injection ports. 
Using current rocket engine technology, a chamber with zirconium- copper, 
channel wall cooled combustor and tubular nozzle regenerative ly cooled with 
hydrogen was designed for a 40-to-l area ratio nozzle. The predicted vacuum 
specific impulse was 1363.3 Ib^sec/lb^. 

An optical train analysis of the 5000-kw configuration resulted in a 90-degree 
orientation on the incoming laser beam and the line of thrust to provide the 
optimum beam receiving flexibility. 
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INTRODUCTION 


With the development of laser energy progressing as rapidly as it is, 
the possibility of beamed energy transmission becomes closer to reality 
with more efficient beam generation. One of the more promising uses of 
laser energy is the NASA-LeRC laser-heated rocket thruster where laser 
energy is beamed to a space vehicle from a gound, air, or space generating 
station. An optical train on the vehicle receives the beam and focuses it 
within the thruster through a solid window. As shown in Fig. 3, the 
absorption of the laser radiation is to occur via inverse bremsstrahlung 
and form a high temperature plasma within the flowing propellant. The 
propellant for this application is hydrogen possibly with s S-micron carbon 
seed. The carbon seed provides a radiation absorption layer to aid in the 
cooling of the chamber wall. This concept also employs an electric arc 
to initiate propellant gas ionization. The attractive features of the 
laser-heated rocket are: 

1. Energy transfer from a gound or space laser energy 
generating station where weight constraints are not 
limiting to a weight-critical space vehicle. 

2 . Potential of 1000 to 2500 lb£sec/lbm specific impulse 
with increase in vehicle payloads. 

A twelve-month program was performed for NASA-LeRC by Rocketdyne as the 
prime contractor and Physical Sciences, Inc., as a subcontractor. A 
space vehicle application using 5000-kw input laser power was conceptually 
evaluated, and a detailed evaluation and fabrication of a 10-kw experi- 
mental thruster was completed. The objectives of this program were to: 

(1) perform design and analysis of a 5000- kw laser heated rocket thruster; 

(2) design a 10-kw configuration through scaling and design analysis; and 

(3) fabricate and deliver 10-kw experimental hardware (thrust chamber 
assembly, thrust stand, and plasma initiation system). 




The program was divided into two technical tasks a:nd one reporting task. 
Task I was to provide design and analysis of the 5000>kw and 10-kw designs. 
The fabrication of the 10-kw experimental hardware was to be performed in 
Task II. 
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CONCLUSIONS 





! 

I I 

The overall conclusions of the analysis tasks performed in this | 

program included: \ 

s 

i 

1. Delivered specific impulses exceeding 1200 Ib^sec/lbm • 

are attainable based on predicitions. \ 

\ ) 

i \ 

2 . Low laser input power results in small chamber sizes 

a. Carbon seed in hydrogen does not significantly 

influence plasma size. i 



3. 10- kw thruster was adequately cooled using water forced 
convection cooling without carbon seeding. 

4. 10- kw thruster zinc selenide window was adequately 
cooled (433 K or 319 F maximum temperature) using 
10-percent of the hydrogen flow as window film coolant. 

5. For the 5000-kw regeneratively-cooled thruster, a 
carbon seeded hydrogen layer is required to aid in 
cooling the chamber wall in the plasma region. 

a. Low seed gas velocity- to-hot gas velocity 
ratios are required. 


6. 10-kw Optical Train 

a. For minimum blur circle size at focus, a single 
transverse mode laser output beam is required. 



7. The unique low cost concept of using an uncooled copper 
chamber with graphite inserts will allow the experi- 
mental evaluation of different internal chamber contours. 

8. The flowing plasma concept evaluated for both the 10- kw 
and 5000- kw configurations combine the convective and 
radiative heat fluxes in the chamber throat. Perhaps 

a stationary plasma concept which would only have the 
convective heat flux in the throat should be evaluated 
for advantages over the flowing plasma concept. 

9. More detail evaluation of the 5000-kw configuration 
should be performed which would include a detailed 
analysis of the carbon seeded hydrogen flow injection, 
plasma initiation, optical train cooling, and the engine 
system (mounting of the thrust chamber and optical train 
to the vehicle) . 


TASK I: ANALYSIS AND DESIGN 

This task consists of six analysis subtasks and two design 'subtasks. Tlie 
analysis subtasks defined the laser-heated rocket chamber configuration, size, 
performance, cooling and optical train requirements for the 10-kw and 5000-kw 
configurations. The two design subtasks provided the design layouts of the 
two configurations and provided detail design drawings for the fabrication 
of the 10-kw experimental hardware. Also, detail design drawings of the opti- 
cal train for the 10-kw thruster were performed as part of the optical train 
analysis. 

PERFORMANCE ANALYSIS 
Theoretical Perfo i yian ce 

Theoretical performance levels which enable the assessment of a rocket 
thruster's performance potential are chemical equilibrium or shifting, and 
frozen values. Chemical equilbriura represents an upper limit in achievable 
performance and is the reference condition used in JANNAF performance metho- 
dology. All applicable thruster performance losses are subtracted from the 
©quilibrium value to determine the delivered value. For chemical equilibrium, 
the propellant is assumed to maintain chemical equilibrium (i.e., chemical 
species continually change) as it progresses through the throat and the nozzle 
expansion. In actuality, the nozzle flow has finite reaction rates aad does 
not occur instantaneous as assumed in the chemical equilibrium value. For 
frozen data, the gas composition is fixed at tliat in the combustion chamber 
and represents a lower performance limit. 

To predict the delivered specific impulse of the 5000-kw (40:1 area ratio) and 
the 10-kw 1-014 X 10^ N/M^ or 14.78 psia) designs, parametric theoretical 

specific impulse data for various chamber gas temperatures, carbon concentrations, 
and cesium concentrations were generated for a chamber pressure of 3.45 X 10 N/M 
or 50 psia. Both equilibrium and frozen data were obtained. Since if carbon is 
used it will be injected as a solid and may or may not chemically react with the 
hydrogen, the influence of non-reacting and reacting carbon was evaluated. 
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For the hydrogen- carbon mixture, the addition of the high molecular weight 
carbon degrades the theoretical equilibrium specific impulse as shown in 
Fig. 4 and 5. This degradation becomes greater at the higher chamber tempera- 
tures due to the greater amount of dissociated hydrogen at these temperatures. 

Also shown in Fig. 4 and 5, the non-reacting or solid carbon state resulted 
in a 70-lb^sec/lbm (approximately 5 percent) lower specific impulse than if 
the carbon is assumed vaporized and reacting at an average chamber temperature 
of S555.6 K (10,000 R) and for the 40:1 area ratio. However, at high tempera- 
I tures (3888.9 K or 7000 R and above), the carbon has been theorized to I 

i vaporize and react with hydrogen and no longer remain as a solid. Therefore, 

I more realistically the maximum specific impulse difference between solid 

I (non- reacting) carbon and reacting carbon will be less than 30 Ib^/sec/lbm 

(2 percent) at 3888.9 K. 

i 

i Comparing Fig. 4 (40:1 area ratio) and Fig. 5 “ 1.014 X 10 N/M or 

j 14.7 psia), the lower area ratio condition resulted in a 30-percent lower 

theoretical specific impulse. The influence of nozzle area ratio is shown in 
Fig. 6. A specific impulse increase of approximately 7 percent can be 

( i 

; achieved with an increase in area ratio from 40 to 200. The carbon con- 

: centration is denoted by C which represents the weight fraction. 

I c 

i i 

For the 5000-kw design to achieve a vacuum specific impulse greater than a j 

I I 

1000-lb^-sec/lbm, the average chamber gas temperature will have to exceed | 

3333.3 K (6000 R) as shown in Fig. 7. The increase in the slope of the | 

specific impulse versus temperature curve above 2777.8 K (5000 R) is the ■ 

i result of a significant increase in dissociated hydrogen at these temperatures 

i (Fig. 7 and 8). For pure hydrogen, this theoretical equilibrium vacuum | 

j specific impulse variation with temperature up to 15,000 K (27,000 R) is 

j presented in Fig. 9. 
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Figure 6. Influence of Area Ratio on Tieoretlcal Equilibrium Vacuum Specific Impulse 
(555S.6*K or IO.OOO*R) 

































THEORETICAL SPECIFIC IHPULSE, LB.-SEC/LB, 


• 1000 




CARBON, FRACTION 

V- ■ -j . ’ ’ ; . 

, ... ... 

. BY WEIGHT 

t A 


hydrogen-carbon 

Z.‘t5 X 10^ N/M^ 

*' (50 PS I A) 

QUI LIBRIUM 

1.014 X 10^ 

N/m2 (I4J_PSIA 




I 


TEMPERATURE, R 


TEMPERATURE, K 


TheorotiCAl EqullibxluB Specific Iqpulse Varietion 
with Chaaber TeaDerature and Carbon Concentration 

^**exit " 1*014 X loS n/M* or 14.7 psia) 






PULSE » L8f-SEC/LB 





The difference between equilibrium and frozen specific impulse is presented 
in Fig. 10 and 11 as a function of chamber temperature and carbon concentration. 
This difference represents the potential reaction kinetic loss which is de- 
pendent on the propellant conditions (chamber pressure, temperature, and 
carbon concentration), chamber size, and nozzle area ratio. The lower area 
ratio (approximately 1.2:1) of the 10- kw design tends to lower kinetic losses. 

However, the small size of this design v/ill increase the kinetic loss. There- 
fore, the frozen specific impulse variation with chamber temperature and car- 
bon concentration was presented in Fig. 12 and 13 for an ambient exit pressure. 
As shown in Fig. 10 and 11, the difference between equilibrium and frozen 
specific impulse decreased with increase in carbon concentration. 

The addition of cesium to the hydrogen-carbon mixture can enhance the plasma 
formation and initiation. The influence of adding cesium to the hydrogen- 
carbon gas at a chamber temperature of 3888.9 K (7000 R) is shown in Fig. 14 
and 15. Similarly, for the carbon the addition of cesium tends to decrease 
the theoretical specific impulse and the difference between equilibrium and 
frozen values. 


Typical equilibrium gas composition distributions with the carbon reacting 
are presented in Fig. 16. The distribution starts at the injector and pro- 
gresses downstream to the nozzle exit (40:1 area, ratio). As shown in Fig. 16, 
the majority of the gas is composed of H and H 2 . The recombination of H to 
11^ in the nozzle is evidenced by the decrease in H and the increase in H 2 . 
Also, the equilibrium condition predicts the formation of carbon at an area 
ratio of 12.5:1. However, it is anticipated that the reaction rate of form- 
ing solid carbon is so slow that solid carbon will not form. 
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CARBON, FRACTION BY WEIGHT 


Figure 12. Theoretical Frozen Specific Impulse Variation with 
Carbon Concentration and Chamber Temperature 
- 1.014 X 105 n/M^ or 14.7 PSIA) 
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Pisure 13. Theoretical Frozen ^ecific Impulse Variation 
with Chtnber Teaperature and Carbon 
Concentration (P . i.oU x lo5 h/¥T or U.7 PSIA) 
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CESIUH. FRACTION BY WEIGHT 


Figure IS. Specific Inqpulse Difference (Equilibrium and Frozen) 
Variation with Cesium and Carbon Concentration 


^‘^chamber * 3888.9 K or 7000* R) 







SPECIES CONCENTRATION 











Delivered Specific Impulse 


P^erformance Loss Mechanism. The performance loss mechanisms associated with 
the laser-heated rocket are similar in principle to those which occur in a 
conventional chemical liquid rocket thrust chamber. Therefore, the procedures 
used in calculating the performance losses can be fashioned after JANNAF 
methodology. 

Propellants (hydrogen, carbon and/or cesium) enter the thrust chamber through 
an injector. The laser beam sustains the plasma vaporizing the carbon and/or 
cesium and heating the remaining gas. The performance loss mechanisms of the 
laser-heated rocket are presented in Fig. 17. Any deviations from complete 
homogeneous mixing and heating to equilibrium can be referred to as energy 
release losses (energy release efficiency, The majority of the plasma 

radiant energy will be absorbed by the carbon-seeded hydrogen; however, some 
radiation will be lost through the window (plasma radiation loss, n ). 

Chemical reactions continue to occur in the subsonic region as well as in the 
supersonic nozzle. In the expansion process (nozzle), deviations from the 
local chemical equilibrium condition are losses associated with reaction 
kinetics (reaction kinetic efficiency, n,J . Also in the physical expansion 
process, flow velocity components normarto the desired direction of thrust 
will occur. The loss in momentum and attendant nonuniform expansion is re- 
ferred as a divergence or geometric loss (geometric efficiency, n ) . A 
boundary layer is formed in the region adjacent to the chamber wall due to 
the presence of the wall. Losses in available thrust caused by momentum and 
heat transfer to the wall are termed boundary layer losses (boundary layer 
efficiency, The existence of a nonuniform temperature or carbon and/or 

cesium concentration distribution at the nozzle exit can result in a perform- 
ance loss (nonuniform temperature and concentration efficiency, n ). This 
latter influence can be determined using one of two methods. One approach 
is to use a stream tube analysis and account for the radial distribution of 
mass flow, gas temperature, carbon concentration, and the different final 
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FIGURE 17. PERFORMANCE LOSS MECHANISMS 


stream tube area ratios, Another more simplified approach is to utilize an 
average overall gas tenqpcrature and essentially one-dimensional flow. 

Nozzle Conto ur Des ign . 

10-KW Chamber. Since the 10-kw chamber was to be tested essentially at 

sea level, the nozzle area ratio was sized to have the nozzle flow full with- 

5 2 

out internal flow separation. Based on the chamber pressure of 3.45 X 10 N/M 
(50 psia) and combustion gas properties, an area ratio of 1.6-to-l resulted. 

Also since any performance gained from a contoured bell nozzle would be small for 
this low area ratio, a simple 15-degree conical nozzle was selected. 

5000- KW Chamber To provide an efficient and rapid analysis of the 

5000-kw chamber, the nozzle contour in normalized coordinates was designed 
by analyzing the nozzle contour over a range of conditions. This was accom- 
plished by selecting a parabolic nozzle contour for the 40:1 area ratio, 80- 
percent length nozzle. This type of contour closely approximates an optimum 
contour in geometry and performance. For propellant combinations with high 
temperatures and significant amounts of hydrogen, reaction kinetic loss can 
be significant. This kinetic loss can be minimized through proper nozzle con- 
touring which can be accomplished simply by increasing the upstream and down- 
stream throat radii of curvature. As shown in Fig. 18, this tends to reduce 
the flow expansion rate which in turn reduces the kinetic loss. 

An 80-percent length, 40-to-l area ratio nozzle contour designed with a 
downstream throat radius of curvature of 2, is presented in Fig. 19. 

The conqplexity of radial varying temperatures and carbon concentration pre- 
cludes the analysis of the supersonic flowfield using full reacting flow. 

Past analyses conducted with chemical rocket thrust chaad>ers indicated that 
a constant gamma analysis using the chemical equilibrium throat gamma (y*) 

Closely approximated the reacting flowfield analysis. The variation of the 
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Pigim 18. In£lu 0 nc 0 of Upstrom Throot Rodlut of Curvoturo on 
Plow Expansion Rate 
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Typical Hozzle Contour 





equilibrium throat gaimna with an average chamber temperature and carbon con- 
ccntr?.cion is presented in Fig. 20 • For a temperature range of 3888. 9K (7000 R) 
to 9555.6 K (10,000 R) the gamma varied from 1.16 to 1.28. 

The nozzle contour shown in Fig. 19 was analyzed for these two gammas. First, 
c transonic solution was performed to describe the flowfield in the 0.8 to 
1.2 Mach number region. This was accoaq;>lished using a computer which approxi- 
mates the parabolic partial differential equations using a 33-term power 
series expansion. The results of this program were fed into the nozzle 
analysis program which performs the characteristic net solution for the re- 
maining supersonic flowfield. Results of this analysis include the geometric 
efficiency and the flowfield characteristic net which could be transfort ;d 
into stream tubes. 

Reaction Kinetic Analysis . To perform a reaction kinetic analysis, the perti- 
nent chemical reactions and their reaction rates must be defined. For the 
hydrogen- carbon propellant combination, two sets were defined. As shown in 
Table 1, one set is a simple set with only three reactions and the other is 
a more complete set with ten reactions. The third body efficiencies for the 
reactions containing a third body (reactions I, 2, and 3) are presented in 
Table 2. The difference in reaction kinetic efficiency for these two re- 
action sets was evaluated for an arbitrarily selected condition of an average 
chamber temperature of 3888. 9K (7000 R) and 20-percent carbon by weight. 

Using the one-dimensional kinetic (OOK) computer program, the difference in 
kinetic efficiency was only 0.2 percent. Therefore, the simple reaction set 
was used. 

10-KW Chamber Delivered Specific Impulse . The 10-kw chamber contour, with a 
1.6:1 area ratio, 15-degree conical nozzle, was analyzed using the Rocketdyne 
Boundary Layer Computer Program which utilizes an integral method of solving 
the momentum and energy equations (Table 3) to determine the boundary 
layer loss. A laminar boundary layer was assumed. The results of this 
computer program were also used in determining the convective portion of 
the chamber heat flux in which the boundary layer recombination influence 
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TABLE 1. CARBON-HYDROGEN REACTION SET 

























10- KW CHAMBER PERFORMANCE 




mmm 


1232.5 lb, sec/lb 



was neglected. The ODK COne-Oimensional Kinetic) analysis resulted in a re- 
action kinetic efficiency of 0.9765 for the 10-kw chamber. The geometric 
efficiency of the 1.6:1 area ratio, 15-degree cone was 0.9662. The simple 
source flow solution of a 15-degree conical nozzle yields a geometric 
efficiency of 0.9830, however, due to influences of throat radius of curvature 
and the low area ratio the geometric efficiency results in a lower actual value. 
These predicted efficiencies were used in the following simplified JANAFF 
equation. 


I 

^del 


I 

*theor 

vac 


"c* "g \ V ■ - "bl’ 


( 1 ) 


where I = Delivered vacuum specific impulse 

®del 

Ig = Theoretical vacuum specific impulse 

theor at injected condition 
vac 

= Characteristic velocity efficiency, or n j,. 
Energy release efficiency 

= Geometric or two-dimensional efficiency 
= Reaction kinetic efficiency 

Ohl ~ Heat loss efficiency (includes and r)^) 

rigL “ Boundary layer efficiency 


For this analysis, an effective combustion gas temperature was used rather 
than the stream tube method; therefore, the non-uniform temperature efficiency 
(n ,. ) is 1.0. A discussion of the method of computing this effective tern- 
I'-erature is presented in the Flow Analysis section. Since an effective combus- 
tion gas temperature was used which accounted for heat losses, the values of 
and were 1.0 and the theoretical vacuum specific impulse of 1465 Ib^ 
sec/lb^ at the effective temperature of 5556 K (10,000 R) was used. The de- 
livered specific impulse of the 10-kw chamber was computed to be 1232.5 Ib^ 

sec/lb as shown in Table 3, which is well above the contract specific 
m 

impulse goal of 1000 Ib^sec/lb^. The predicted delivered thrust is 0.387 
Newtons (0.087 Ibf). 
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5000- KW Chamber Delivered Specific Impulse. Using the same procedure out- 
lined for the 10-kw chamber, the boundary layer and reaction kinetic effi- 
ciencies were obtained as shown in Table 4. The geometric loss was 
determined by performing an aerodynamic analysis of the 40-to-l area ratio 
bell nozzle. Using the computed effective combustion gas temperature 
(Teff = 5556 K or 10,000 R) the predicted delivered specific impulse for 
the 5000-kw chamber was determined to be 1363.3 Ib^sec/lb^. The predicted 
delivered vacuum thrust is 351.1 Newtons (7 9.2 Ib^) . 




PLASMA ANALYSIS 


TTie detailed plasma analysis for this program were performed by Physical 

Sciences, Inc. (PSI) and presented in Ref. 1 . The pertinent results are 
presented in this section. 

Plasma Absorption 

Using the hydrogen equilibrium species composition predictions of Patch over 
the pressure range of 1 to 10^ atmospheres and temperature range of 6000 K 
(10,800 R) to 18,000 K (32,400 R) , the resulting absorption coefficient for 
10.6 pm radiation (normalized by pressure) is shown in Fig. 21 . Since 
neutral and ionic inverse Bremsstrahlung exhibits different pressure depend- 
ence, the individual absorption coefficient, rather than the sum, has been 
plotted in order to facilitate interpolation between curves. Tlierefore, the 
total absorption coefficient is the sum of the two curves. For example, if 
an absorption coefficient of 1 cm'^ (2.54 in.'^) was desired at a pressure of 
three atmospheres, an operating temperature of approximately 13,000 K 
(23,400 R) would be required. 

Observing Fig. n , the absorption coefficients are very small at temperatures 
less than 10,000 K (18,000 R) which represent long plasma lengths. To oper- 
ate at lower temperatures, one possibility is to seed the hydrogen with a 
readily ionizable gas such as cesium. Data were obtained for several cesium/ 
hydrogen mixtures; a typical variation with temperature and pressure is 
illustrated in Fig. 22 which exhibited an interesting trend with temperature. 
The initial rise in absorption coefficient at 3000-4000 K (5400-7200 R) is due 
to incomplete cesium ionization. The cesium is typically fully ionized at 
temperature greater than 4800 K (8640 R) . The second increase occurring at 
temperatures of approximately 8000 K (14,400 R) is due to the onset of sig- 
nificant ionization of the hydrogen. Tlie absorption coefficient may be in- 
creased by increasing the cesium concentration; however, this can only be 
done at the cost of a severe weight penalty since the cesium/hydrogen atomic 
mass ratio is approximately 133. 
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INVEfSE BREMSSTRAHLUNG ABSORPTION COEFFICIENT FOR 
HYDROGEN PLASMA Vs. TEMPERATURE AND PRESSURE, x -10.6 uM 













One concept is to provide a layer of seeded gas containing carbon particles 
as a plasma radiation absorption flow to provide a means of preventing the 
plasma radiation from reaching the wall. Since the coolant and plasma flows 
will mix to some degree, carbon particles may be along the laser path and 
could provide an absorption mechanism for the laser radiation. 

The absorption coefficient for small carbon particles with radii within one 
or two orders of magnitude of the wavelength of the incident radiation was 
determined through use of Mie theory. The required parameters for evaluation 
of the Mie equations are the complex index of refraction of the material at 
the wavelength of interest along with the ratio of particle radius to wave- 
length. Unfortunately, the complex index of refraction of carbon at 10.6pm 
is not well defined, Mie theory prediction were performed for three repre- 
sentative values of the complex index of refractions; (a) m ■ 4.0 - i 0.3, 
representative of measurements made of soot both at room temperature and at 
2200 K; (b) m “ 4.0 - i 1.0, representative of room temperature measurements 
of pyrolitic graphite, where the imaginary part is uncertain to within a 
factor of two; and (c) m = 4.4 - i 4.0, corresponding to a theoretical predic 
tion for graphite particles at 2100 K. While it is not expected that the 
complex part of the index of refraction of graphite would increase by a 
factor of 4 between 300 K and 2100 K, this last case (c) has been included to 
demonstrate the possible temperature variation in the absorption coefficient. 

Since the interest of this work is in defining the absorption length for a 
given particulate mass loading, it is convenient to express particle absorp- 
tion properties in terms of a, the volume absorption coefficient in cm"^. 

The absorption cross section for a particle Is defined by 

Op • a Vp (2) 

where Vp is particle volume and the total absorption length for a given 
particle density, Np, is 
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or in terms of particle mass loading/cc, m 


L = (om /p )'^ (4) 

P P 

shere is the particle specific density. Tims, o may be used to directly 
I elate absorption length to mass loading. Note that the most efficient 
absorption/unit mass occurs at the largest value of a. 

Predicted values of a for the three aforementioned values of the complex 

index of refraction are shown in Fig. 23, for particle radii between lO”^ - 
+2 

10 pm. It can be seen in all three cases that the most efficient absorption 
occurs for particle sizes of 1 - 2 um. Tlie decrease in a with increasing 
particle size is due to optical thickness effects. 

lor purposes of example, if a is taken to be 6 x 10^ cro‘\ then from Eq. 4 a 
mass loading of 3.7 x 10 gm/cc would result in an absorption length of 1 cm, 
wnere p^ has been taken to be 2.25 gin/cc. (For comparison, 3 atm. of H 2 at 
room temperature would have a density of 2.4 x 10*^ gm/cc.) It should be 
pointed out that th*. absorption efficiency of the graphite particles would 
<lec.';ise at higher temperatures because of vaporization effects since the 
boiling point of graphite is 4500 K. The magnitude of this effect would, of 
course, be dependent on flow conditions. 

Plasma Core Sizing 

lUscussions with PSl indicated that the most realistic plasma core concept was 
the continuous flowing plasma flow pattern (Fig. 24). Therefore, the effort 
Concentrated on this approach. 









Figure 24. Laser-Heated Rocket Thruster — Continuous Flow Plasma Concept 



The plasma core analysis utilized a one-dimensional laser supported combustion 
(LSC) wave analysis and an available nozzle computer code to define the plasma 
core properties. 

An LSC wave analysis analogous to that of Jackson and Nielsen (Ref. 2) was 
performed to estimate the propagation velocity of the plasma as part of 
Contract NAS3-19695, Laser Heated Rocket Absorption Study. In order for the 
high temperature plasma to remain stationary within the chamber, the front of 
the "heated zone" (plasma) must propagate at a velocity equal to but opposite 
to inlet flow velocity (Fig. 24). The variation of the hydrogen plasma pro- 
pagation velocity ("cold gas") with laser flux or intensity is presented in 
Fig. 25. In addition, theoretical results (Ref. 2) and experimental data 
for air are presented for comparison. 

The variation of the initial velocity forward edge of ("hot gas") plasma 
with laser intensity is shown in Fig. 26 for a chamber pressure of 3 atm 
(~50 psia) and temperatures of 10,000 K (18,000 R) and 12,000 K (21,600 R) . 
From the nozzle code the variation in inlet velocity with throat area-to- 
inlet area ratio is defined. Therefore, for a given laser intensity the 
intersection of the LSC wave theory curve and the nozzle code curve (Fig. 27 ) 
defines the throat area-to-inlet area ratio. Performing this for laser in- 
tensities from 10^ to 10^ watts/cm^, the ratio of "co'e" throat area-to-inlet 
area with laser intensity shown in Fig, 2 8 was determined. As shown in this 
figure, the area ratio (inverse of contraction ratio) is not significantly 
influenced by the plasma temperature. The laser intensity threshold to main- 
tain a sustained H, plasma is approximately 10^ watts/cm^. As the laser in- 
tensity is increased to greater than 10^ watts/cra^ laser supported detonation 
wave would form. For the 10-kw configuration the variation of the area ratio 
with plasma inlet area is obtained from the data presented in Fig. 28, as 
shown in Fig. 29. The plasma inlet diameter varies from approximately 0.1 cm 
(0.0394 inch) to 0.4 cm (0.158 inch) depending on the laser i.'‘tensity. The 
corresponding plasma throat diameters vary from approximately 0.08 cm 
(0.0315 inch) to 0.1 cm (0.0394 inch). 
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Plasma Inlet Velocity Variation With Laser Intensity 
and Plasma Temperature 
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Plasma Throat Area-to-Inlet Area Ratio Variation with Laser 
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Selecting a laror intensity of 3*7 x 10 watts/cm , the flow analysis for 
the 10-kw configuration resulted in the pressure and absorption coefficient 

X 

profiles shown in Fig* in x space where r =^/ dx. The plasma static 
temperature profile in real space is shown in Fig. 31. The actual plasma 
core dimensions are presented in Fig. 32. Although this does not include 
the outer buffer gas flow, the relatively small dimensions of the plasma 
core will necessitate a funnel-type chamber configuration to provide suffi- 
cient dimensions in the window and injector region. 

In the process of sizing theVlO-kw plasma core, PSI determined an approxi- 
mate scaling law of plasma diameter with laser beam power. 

Diameter a ^PoweT (5) 

Using this law, the plasma core size for the 5000- kw configuration was de- 
termined (Fig. 33). The diameter of the plasma varied from 4.13 cm (1.626 in.) 
at the plasma inlet to 1.98 cm (0.780-in.) at the geometric throat. 

Plasma Radiation 

The high temperature plasma contained in the chamber loses energy through 
radiative transfer to the chamber walls. This power loss or wall heat flux 
contribution was estimated by PSI and was used to determine wall cooling re- 
requirements. The plasma radiation loss from the plasma for the 10-kw hardware 
was computed to be an average of 1200 watts/cm (7.34 BTU/in, sec) radiation 
heat flux and is primarily from the ultraviolet wavelengths. For the 5000- kw 
configuration, the plasma radiation loss is from both the ultraviolet (UV) 

and infrared (IR) wavelengths which increased the radiation heat flux to 
2 2 

6400 watts/cm (39.13 BTU/in. sec). 

Plasma Radiation Absorption by Particles 

Small carbon particles, smaller than 1 ym, are very efficient absorbers of 
radiation in the visible and infrared wavelength regions, approximately acting 
as black-body absorbers. A large fraction of the hydrogen plasma radiation 
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X, cm (FROM PLASMA INLET) 


Figure 31. Static Temperature Distribution for 10-kw Power Input 
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Figure 32. 10-kw Thruster Plasma Core Dimensions 
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will unfortunately originate in the far UV where little is known about the 
absoiptivc properties of carbon. In the following analysis it was assumed 
fhat the carbon particles absorb as black bodies at all wavelengths. 

'"he scenario considered is a cylindrical plasma of radius r^ surrounded by 
a particle seeded coolant flow of thickness Ar. The particles are taken to 
i>e spherical with radius and to liave a constant number density Np in the 
'-’oolant flow. Each particle has an absorption cross section of itR ^ and 

P 

thus the optical depth, t, across the coolant layer is defined by 

f A r (6) 

■\n optical depth oi 3 is sufficient to absorb 95 percent of the plasma radia- 
tion and this condition is met when the particle number density is given by 


•p 2 2 

^ TT Rp A r Rp Ar 

It is perhaps more interesting to examine this quantity in terms of the total 
mass/cc of particles in the coolant flow, M^, 

"p ■ ■■'p '"p^ (8) 


where is the specific density of carbon which is ~ 2 gm/cc. Combining 
Hq. 7 and 8 results in 


M 




(9) 


Equation 9 can readily be extended to deduce the mass flow of particles 
required to absorb the plasma radiation. 





p c coolant 


( 10 ) 


.'15 


where 


Coolant * Tj Ar A Tj ) 

and is the coolant flow velocity. Since it xs desirable to maintain a con 
figuration such that Ar/r^ « l, the approximation 

Coolant" ^ ’^l (12) 

is made. Combining Eq. 9, 10, and 12 results in 


Mp = 8 Pp Uc % ^1 “ 48 Rp rj 


(13) 


Thus it can be seen from Eq. 13 that as lonq as Ar « r^. the particle mass 
flow required to absorb the plasma radiation is independent of the coolant 
flow thickness and linearly proportional to both the particle and plasma 
radii. Since in the interest of optimizing thrust one wishes to minimize the 
coolant flow, it is clear that the smallest particle size, consistent with 
tlie assumption of blackbody absorption, should be used. 


Another interesting feature is that since the required coolant flow is linearly 
proportional to plasma radius, the ratio of coolant flow to plasma flow is 
size dependent, i.e.. 


M 


8 Pn^nU, 

E-E.L 


M 


plasma ^plasma *^1 *^pl 


(14) 


where Upj is the plasma flow velocity. For example, for a 14,000 K hydrogen 
plasma at a pressure of 3 atm 


M 


10 


M. 


plasma 


R U 

P c 

r,U , 
1 pi 


(IS) 


Note that relationship (15) imples that for small devices, i.e., 1 cm, 

the coolant mass loading would of necessity be orders of magnitude larger 

than that of the plasma for reasonable particle sizes unles U « U . 

c pi 

For the 5000-kw chamber assuming a 1 pm carbon particle and a 1.7 cm plasma 
radius, the required velocity ratio was approximately 1:140. This may or may 
net represent a realistic situation for a single point seeded carbon injection 
but could be for a multi- injection scheme. 

Plasma Arc Initiatio n 

Three power supply systems shown in Fig. 34, were evaluated for potential use 
with the 10-kw cxpei’imental hardware to initiate the gas breakdown. 1'esJa 
coil initiation has proven successful in earlier work initiating breakdown 
in a hydrogen atmosphere when laser power (pulsed) is approximately 10^ watt/ 
cm . However, this device may still not produce the high electron densities 
required of lO^Vcc. A power supply like the self-firing spark gap shown in 
Fig. 34, is capable of producing these electron densities and delivering 2000 
amps peal current. As shown in Fig. 35, (Ref. 3), a peak current of 2000 amps 
produced the required 10 /cc electron density in hydrogen. The self-firing 
spark gap device would arc whenever the capacitor discharged and the exact 
time may not be predictable. 

I he most versatile and powerful system would be the series injection triggered 
spark gaj . This s> stem would provide a range of peak current of 1000 to 5000 
;4jiips with the 10-ku power supply. Ihe capacitor is charged and the series 
injection trigger i.iodulc would cause the capacitor to dump at the exact speci- 
tied tirae. This would permit the exact time sequencing of the laser bea»», 
|iti»llcrt, and initiating arc. Iherefore, this system was selected as shown 
in I’ig. 
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Candidate Power Systems 







SELECTED 10-KW CHAMBER PLASMA ARC INITIATING POWER SYSTEM 
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FLOW ANALYSIS 


10-kw ChaR'.ber Sizing 

The procedure used in sizing the 10-kw and 5000-kw chamber is schematically 
shown in Fig. 37 through 39. The sizing included potential flow (not in- 
cluding boundary layer thickness) sizing of the plasma core and chamber 
throat (Fig. 37). The latter was sized assuming an average gas temperature 
and carbon concentration. From laser supported combustion (LSC) wave theory, 
the flow velocity at the plasma Inlet was specified which defined the chamber 
diameter at this point. Then a heat balance between the plasma, the colder 
buffer gas, and the chamber wall cooling was performed. The schematic of 
the heat balance performed is presented in Fig. 38 and 39. The heat balance 
consists of boundary layer analyses for different freestream conditions and 
a conduction analysis of the flow between the plasma and the chamber wall. 
This heat balance results in heat fluxes, freestream conditions, and an 
average gas temperature. Now accounting for the heat lost to the chamber 
vvall and through the window, the effective average gas temperature is com- 
puted. This gas temperature is iterated until the assumed average gas 
temperature is close to the calculated value. 

The total propellant flowrate may be computed by 


, (P) (0.948) 

‘''total ♦ I*’/ sec 

A H 


(16) 


and the potential or inviscid flow (i.e., not including boundary layer 
influences) chamber throat radius is given by. 





0.948 P c* 

theor 

- i 


IT g AH P^ 


inch 


(17) 
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Heated Rocket Chamber Sizing 







^ 9 * StreaiPtubf Heat Balanc 







where P 

A ii 
0.948 
c* 

g 

P 

c 


Laser beam power input, kilowatts 
Propellant enthalpy rise, Btu/lb 
Conversion from kilowatts to Btu/sec 
Theoretical characteristic velocity, ft/sec 
Gravitational constant (32.174 ft/sec^) 
Chamber pressure, psia 


The variations of these two parameters with average gas temperature and carbon 
concentration for the 10-kw configuration are presented in Fig. 40and 41. 

Pure hydrogen at an average gas temperature of 5555.6 K (10,000 R) resulted 
in a 0.0541 cm (0.0213 inch) inviscid flow throat radius and a total propel- 
lant flowrate of 0.0000704 Ib/sec. For this condition, the plasma core flow 
would be 31.5 percent of the total flow. 


Boundary layer analyses were performed using the Rocketdyne Boundary Layer 
Analysis Computer program assuming a 10 -kw chamber sized for an average gas 
temperature of 5555.6 K (10,000 R). As shown in Fig. 42. the resulting con- 
vective heat flux (maximum and average) on the chamber wall with variation in 
the boundary layer edge temperature was obtained. An eleven nodal one- 
dimensional them.al conduction model of the gas was constructed and included 
the variation in the g&s thermal conductivity (Fig. 43). Using an average 
buffer gas (hydrogen) thickness for the chamber sized for 5555.6 K (10,000 R) , 
the conduction analysis was performed for various boundary layer edge tempera- 
tures. The resulting conductive heat flux variation with boundary layer edge 
temperature is shown in Fig. 44. In addition the resulting temperature pro- 
files through this gas layer are presented in Fig. 45 . 

The heat flux absorbed by the buffer gas was determined using two methods of 
computation: 

1. From boundary layer and conduction analysis 

2 , From propellant properties and conduction analysis 
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Figure AO. Hot-Gas Flowrate Variation with Gas Temperature and Carbon 
Concentration (10 KW Configuration) 
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Figure 41. Chamber Throat Radius Variation With Gas Temperature and 
Carbon Concentration (10 kw Chamber) 
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Figure 45. Temperature Profiles from Conduction Analysis 
(10 KW Chamber) 
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The conditions were varied until the heat flux absorbed computed from the two 
calculational methods were equal as illustrated in Fig. 46. The pure hydrogen 
was assumed to be transparent to the plasma radiation. As shown in Fig. 4o, 
the heat balance resulted in a 4222.2 *'f7600 R) boundary layer edge temperature. 
The average overall temperature of 7972.8 K (14,351 R) was calculated for the 
no heat loss condition. The average convective heat flux plus the radiation 
heat flux to the cliainber wall (9.75 Btu/in^ -sec or 1591.4 watts/cra j repre- 
sents approximately a 20-percent heat loss. Assuming a 5- percent heat loss 
through tlie window, the effective average overall gas temperature was computed 
to be 5888.9 K (10,600 R) which is amazingly close to the 5555.6 K (10,000 R) 
assumed. 

The dimensions of the resulting 10-kw chamber are shown in Table 5 and Fig. 4". 
Considering the flow heat balance ;.nd the boundary layer influences, a geo- 
metric throat diameter of 0.1168 cm (0.046 inch) was determined. This con- 
sists of 0.1082 cm (0.0426 inch) for the potential flow and 0.00864 cm 
P.0034 inch) for the boundary layer thickness. The geometry of the region 
n''ar the window and injector illustrated in 1 ig. 47was approximated based on 
window requirements, flow recirculation tendencies, and optical interface 
difficulties . 

10-kw Combustion Chamber (Subsonic) Flow Fat tern 

The subsonic flow pattern of the 10-kw chamber was estimated by observing past 
water test results and flow pattern interpretation of actual hot-firing rocket 
chambers. In subsonic flow with fixed fluid injection locations, the flow 
recirculation patteni will be independent of tlie injection velocities. An 
increase in injection velocities will merely increase the vortex velocities. 

The predicted subsonic flow pattern for the 10-kw chamber is presented in 
Fig. 48. Basically two vortices will exist. One will be located near the 
window flowing counter-clockwise and the other located near the main propel- 
lant injection flowing clockwise. 
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I'O provi.le a defined flow of propellant along the wall of the chamber, the flow 
will he injected parallel to the chamber wall at a velocity the same order of 
magnitude as the equivalent injector annulus velocity. This injector annulus 
velocity IS the velocity of the injected propellant through the annulus formed 
bv the chamber wall and the inner edge of the injector lip (Fig. 48). A main 
propellant injection velocity of 152 cm/sec (5 ft/sec) and a window coolant 
velocity of 40 cm/sec (1,5 ft/sec) were selected. 


10-kw Chamber/Facility Flow Stability 

A phenomenon which may occur in a test chamber is a coupled chamber-to-test 
facilit} flow instability caused by a low chamber injector pressure drop. A flow 
disturbance occurring within the chamber propagates into the facil-’ty line which 
decreases the flowrate into the chamber. The chamber pressure then decays and 
the flowrate Increases resulting in an unstable flow condition. Past experience 
has indicated that an injector pressure drop of approximately 10 percent of 
chamber pressure eliminates the unstable condition. Due to the small flowrate of 
the 10-kw chamber, the only feasible method of creating this magnitude of hy- 
draulic resistance was the use of an orifice consisting of a plate of rigimesh. 

50t)0-kw C hamber Sizing 

The sizing of the 5000-kw chamber was accomplished using the same procedure 
used for the 10-kw chamber. The variation of the 5000-kw chamber flowrate and 
the potential flow (i.e., not including boundary influences) throat radius with 
average gas temperature and carbon concentration are presented in I*ig. 49 and 
50. Due to the high plasma radiation heat flux (Plasma Analysis Section) for 
this configuration of 6400 watts/cm“ (59.13 Rtu/in“-sec) . a high concentration 
of carbon was assumed of 0..5. It was felt that SO percent carbon represented 
perhaps an upper limit in efficient mixing with injected hydrogen. The heat 

flux results of the boundary layer and conduction analyses arc presented in 
I'ig. 51 and 52. 
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Figure' 50. Chamber Throat Radius Variation with Gas Temperature and 
Carbon Concentration (5000 kW Chamber) 
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Figure 52. Average Conductive Heat Flux Variation with Boundary Layer 
Edge Temperature (5000 KW Chamber) 
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AVERAGE CONDUCTIVE HEAT FLUX, WATS/CM 


An attempt was made at cooling tlie convective lieat flux (approximately 491 

2.2 

watts/cm or 3 Btu/in*'-sec and 20 i)crcent of the radiation heat flux approxi- 

2 . 2 

mately 1309 watts/cm" or 8 Btu/in"'-scc) . This assumed that 80 percent of the 
radiation was absorbed by the seeded carbon layer. However, as will be dis- 
cussed in the Cooling Analysis section the available hydrogen flow was not 
sufficient to cool the chamber. Increasing the amount of radiation absorbed 
by the carbon layer to 90-]>ercent, the chamber could be satisfactorily cooled. 

Performing the heat balance to size the chamber, a 2778 K (5000 R) boundary 

layer edge temperature resulted. With the 10-perccnt radiation heat flux 

value of 654 watts/cm (4 Btu/in"-sec) and an average convective value of 
2 2 

393 watts/cm (2.4 Btu/in -sec), the effective average overall gas temperature 
was computed to be 5555.6 K (10,000 R) . 

The dimensions of the resulting 5000- kw chamber are shown in Table 6 and 
Fig. 53. Considering the flow heat balance and the boundary layer influences, 
a geometric throat diameter of 2.672 cm (1.052 in) was determined. The 
geometry of the region of the window and injector illustrated in Fig. 53, was 
approximated based on laser beam clearance. To provide a stable plasma the 
LSC theory required that a laser intensity nf 3.7 x 10' watts/cm (specifies 
beam diameter at X = -3 cm or -1.18 in.) that the chamber diameter at this 
location be 6.528 cm (2.57 in.), ihcrefore to provide that the incoming laser 
beam clear the thrust chamber wall, the 45.7() degree included angle W'as neces- 
sary. 

Carbon Seeder 


During the initial portion of the program, the carbon was thought to be 
required for the 10-kw experimental hardware to result in a reasonable size 
plasma and to provide cooling of the chamber, therefore an analysis toward 
the definition of a carbon seeder was performed. The carbon was to consist 
of submicron solid particles and was to be mixed with the hydrogen and intro- 


Table 6. 5000-kw Chamber Dimensions 
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TABLE 7. POWDER FEED UNITS 
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Figure 54. Schematic of Candidate Gas-Solid Mixers 


The conceptual carbon seeder system analysis considered the carbon power 

feeder in a high pressure (approximately 3.45 X 10^/M^ or 500 psia) enclosure 

5 2 

and a low pressure (approximately 3.45 X 10 N/M or 50 psia) enclosure. The 
difficulty and cost in design and fabrication of a high pressure enclosure led 
to the selection of the low pressure carbon powder feeder enclosure. As 
shown in Fig. S5> the selected seeder system has the feeder enclosed and 
purged with gaseous nitrogen to eliminate the hydrogen fire and/or explosion 
hazard. The low pressure mixture of hydrogen and carbon is aspirated by a 
high pressure ejector which will tend to further deagglomerate the carbon. 

This mixture represents the carbon seeded hydrogen flow. 

Since submicron carbon particles tend to cake, the carbon powder must be 
thoroughly dried under high temperature and vacuum before being sealed 
inside the gas-solid mixers. 














COOLING ANALYSIS 


10-KW Cooled Chamber 


Forced convection water-cooling of a hot-firing thrust chamber provides a reli- 
able, durable, and continuously operable experimental device. This cooling 

method has been successfully demonstrated on rocket thrust chambers having heat 

2 2 

fluxes as high as 16,357 watts/ cm or 100 Btu/in. -sec which is a factor of 
6 higher than the predicted peak value of the 10-KN chamber. In addition, this 
concept allows the determination of an axial heat flux distribution. Therefore, 
this method of cooling was selected for the 10-KW chamber. 


As shown in Fig. 47> the region adjacent to the plasma core will be water 

cooled using circumferentially machined coolant channels. For the 10-KW 

2 2 

chamber the predicted maximum heat flux is 2500 watts/cm (15.29 Btu/in -sec ). 

2 2 

The predicted plasma radiation heat flux is 1200 watts/cm (7.34 Btu/in -sec ) 
2 2 

and 1300 watts/cm (7.95 Btu/in -sec ) is the convective portion. 

The approximate bum-out heat flux using water as a coolant is presented in 
Fig. 56 along with the predicted maximum heat flux for the 10-KW chamber. 

As shown in Fig. 56» a water velocity of 15.24 M/sec (50 ft/sec ) provides a 
1.7 safety factor on burn-out heat flux and increasing the velocity to 
30.48 M/sec (100 ft/sec ) resulted in a 2.3 safety factor. 


The following water coolant passages were selected: 


Channel width 
Land width 
Channel depth 
Wall thickness 


0.1524 cm (0.060 inch) 
0.3175 cm (0.125 inch) 
0.1524 cm (0.060 inch) 
0.254 cm (0.10 inch) 


The wide land width which is satisfactory from heat transfer standpoint will 
allow individual coolant flow exit lines of reasonable diameters. The indi- 
vidual water exit lines permit the measurement of the individual flowrates and 
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outlet temperatures which allows the determination of the chamber heat flux 
profile. 

With the selected coolant channel configuration, the resulting two-dimensional ly 
calculated wall temperature distributions for a water velocity of 15.24 and 
30.48 M/sec (50 and 100 ft/ sec.) respectively, are presented in Fig. 57 and 58 • 
The computed two-dimensional temperatures outline the cross-sectional of a 
half channel. The maximum wall temperatures of 587 C (1089 F) and 462 C 
(862 F) are satisfactory for operation of the copper chamber. 

A brief analysis of the water coolant pressure drop and flow requironents re- 

c 2 

suited in a 2.965 X 10 N/M (43-psi) pressure drop and 0.907 Kg/sec (2 Ib/sec) 
for 15.24 M/sec (50 ft/sec) and 1.136 X loS/M^ (172 psi) pressure drop and 
1.814 Kg/sec (4 Ib/sec) for 30.49 M/sec (100 ft/sec). 

As shown in Fig. 59, the heat going into the electrode insulator (alumina) 
can only be cooled by the adjacent water passages and any leakage of pressuriz- 
ing gas (hydrogen) along the electrode. A two-dimensional model of the insula- 
tor and the adjacent water passage was constructed. As shown in Fig. 59, 
without gas leakage past the electrode, the insulator would melt ( > 2172 K or 
3450 F) based on the predicted heat flux value. However, a minute leakage of 
3.17 X 10*^ Kg/sec or 0.000007 Ib/sec (l-percent of the main propellant flow) 
cools the insulator to a maximum temperature (Fig. 60) well below the melting 
temperature of the alumina insulator. The alumina was replaced later with 
lavite, a lower cost and more easily machined material. 

Observing Fig. 47, the risks in hot-firing such a chamber contour (required 
from LSC theory and throat area calculations) without previous experimental 
data are apparent. If the resulting plasma is larger than predicted, the 
chamber wall would be exposed to the 12,000 K (21,600 R) temperature plasma. 
Also, if the laser beam moves significantly relative to the chamber, the 
laser beam would hit the chamber wall. These two occurrences could result in 
hardware damage. Therefore, an alternate 10-KW chamber concept was conceived. 
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Figure 59. 10>kw Chamber Arc Initiating Electrode Insulator Well 
Temperature Distribution (Water Cooling Only) 
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10-kfc Uncooled Chamber 


To provide the test program with a higher assurance of success with the 
10- kw experimental hardware^ an additional chamber (a heat sink chamber) 
was designed and fabricated. The additional chamber is a copper he<*t sink 
chamber with graphite inserts (Fig. 61) which provides a highly '/lexibic^ 
short duration test device. This chamber allows for a larger-than-predic cd 
plasma size and is significantly less critical with regard to the r<>lative 
movement of the laser beam with respect to the chamber. Also with the incor- 
poration of the graphite inserts, the internal gas diamoter may be varied to 
evaluate the plasma-buffer gas flow interaction. Also, the water-cooled 
chamber has the flexible capability of being machined to increase the hot-gas 
diameter by a factor of two. The water-cooled chamber provides a test device 
capable of continuous operation and allows measuring of chamber heat flux. 

Analyses of the uncooled chamber included the chamber configurations with 
and without a graphite insert. Also, the feasibility of a plasma viewing 
window was evaluated. 

Without Graphite Insert . A thermal model of a thick wall cylinder was con- 
structed to simulate the uncooled chamber cross section. An effective convec- 
tive heat transfer coefficient was obtained to represent the combined radiation 
and convective heat flux. The copper heat sink chamber was evaluated near the 
plasma inlet (X ■ 2.5 cm or -1.0 inch) and at the geometric throat (X ■ 0.0 cm). 
Using a conservative maximum wall temperature limit of 649 C (1200 F) for the 
copper, the location at X ■ -2.54 cm (-1.0 inch) resulted in maximum test 
duration of approximately 8 seconds as shown in Fig. 62. 

The thermal analysis of the geometric throat was more difficult in that ex- 
cessive CMputer tine would have been required to obtain the tine to a temper- 
ature of 649 C (1200 F). This was due to the small value of the inner radius 
(0.0594 cm or 0.0234 inch). Therefore, data for larger radii were run 
(Fig. 63) and cross-plotted to obtain the 42. S seconds required to reach the 
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649 C (1200 F) . The data for the various radii clearly indicated that as the 
inner radii were decreased (surface area exposed to the hot gas is decreased), 
the time to reach a fixed temperature increased. 

Therefore, without the plasma viewing window, the hot-fire test duration 
limit of the 10-kw uncooled chamber without the graphite insert is approxi- 
mately 8 seconds (Fig. 64). The limiting location was the plasma inlet region. 

Plasma Viewing Window . To enable plasma viewing, photography, and any 
other diagnostics, a transparent viewing window near the plasma inlet could 
be provided as shown in Fig. 65. Two candidate materials (calcium fluoride 
and fused quartz) were thermally evali ’ted. The primary requirement for the 
window was a high transmittance in the visual wave length. Other pertinent 
property data are presented in Table 8. As shown in Table 8, calcium fluoride 
has five times higher thermal conductivity than fused quartz. However, the 
fused quartz is capable of extreme thermal shock and the calcium fluoride has 
poor thermal shock characteristics. 

Tlie results of the thermal analysis of the calcium fluoride plasma viewing 
window are presented in Fig. 66. The analysis indicated that an uncooled 
calcium fluoride window could last approximately 2 seconds and c jjled to 
steady-state maximum wall temperatures of approximately 649 C (1200 F) . How- 
ever, in both uncooled and cooled configurations a thermal gradient exceeding 
204 C (400 F) existed; therefore, due to the poor thermal shock characteristics 
of calcium fluoride and the relatively high thermal gradients which would 
exist, this material was eliminated. 

Similar analyses for fused quartz were performed and, as shown in Fig. 67* the 
maximum test duration was limited to 1 second with or without cooling. The 
extremely low thermal conductivity of fused quartz did not permit appreciable 
heat conduction to provide a cooler wall. However, the excellent thermal 
shock characteristics of fused quartz would most likely allow the 871 C 
(1600 thermal gradier^ predicted. 
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10-KW UNCOOLED CHAMBER THERMAL ANALYSIS 
(COPPER ONLY) 



• TEST DURATION ~ 1 SECOND (LOW THERMAL CONDUCTIVITY WINDOW) 

• WINDOW LIMITS TEST DURATION 



Figure 65. 10 KW Uncooled Chambe 
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Figure 66. 10 KW Chamber Plasma Viewing Window Temperature 
History - Calcium Fluoride 
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Figure 67. 10 KW Chamber Plasma Viewing Window Temperature 
History - Fused Quartz 




Tlierefore, combining the heat sink and plasma viewing window temperature limits 
(Fig. 64), the viewing window would limit the test duration to approximately 
1 second. Therefore, the thermal analysis indicated that a fused quartz view- 
ing window (uncooled) is feasible but severely limits the test duration of the 
uncooled chamber. The value of the information that can be obtained through 
this window in the 1-second duration must be weighed against the increased com- 
plexity incorporating the window and the difficulty in sealing the hot gas 
using the high temperature seals required. 

With Graphite Insert. IVhen the graphite insert is to be used, the plasma view- 
ing window could be replaced with a solid copper plug. Since the graphite ends 
slightly upstream of the geometric throat, the previous analysis of the 
geometric throat location applies with or without the graphite insert. At the 
plasma inlet, the thermal results (Fig. 68) show that the copper portion would 
reach the 649 C (1200 F) limit before the graphite would reach a limit of 
1649C (3000 F) . A composite of the locations analyzed (Fig. 69) indicated 
that geometric throat location would limit the test duration to ^proximately 
42.5 seconds. 

Thermocouples were placed on the uncooled chamber within 1 inch of the copper 
hot gas wall (see design drawings). These will allow the indirect monitoring 
of the gas-side wall tenperature and will enable the detection of an over- 
temperature condition. 

5000- KW Chamber 


A single uppass and a split-flow cooling circuit (Fig. 70) were evaluated for 
the 40-to-l area ratio SOOO-kw chamber. The first circuit is a single i^pass 
circuit with the coolant inlet at the nozzle exit (40-to-l p. *ea ratio). This 
circuit is perhaps the simplest circuit. The entire thrust chamber could be 
fabricated using channels or using the combined construction technique of 
channels and tubes. The second coolant circuit which is slightly more coupli- 
cated would tend to minimize the coolant pressure drop. The same fabrica- 
tion techniques could be applied for this circuit. 



Uncooled Chamber -* Chamber Temperature History with Graphite Insert 




FIGURE 69. 10-KM UNCOOLED CHAflBER THERNAL ANALYSIS 
(COPPER CHAMBER WITH GRAPHITE INSERT) 




Figure 70. ChaabeT/Nozzle Cooling Circuits 
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Preliminary analysis of tube versus channels for the supersonic ponion of the 
thruster (nozzle) from an area ratio of 10 to that of 40 revealed that for 
comparable wall teo^eratures the tubular design resulted in a factor of 15 to 
2o in weight savings. Therefore, due to this obvious weight advantage, a 
tubular nozzle configuration was selected. 

In the evaluation of the coolant circuit, the high coolant bulk ten^erature 
rise which primarily occurs in the plasma containing region of the chamber 
(X»0 to X» -3 cm or X»0 to X- -1.18 in.), reduced combustor coolant flow of 
the split-flow coolant circuit resulted in higher wall tenperatures than the 
single iqppass coolant circuit. Therefore, the single uppass circuit was 
selected. 

An evaluation of how much of the plasma radiation heat flux the regenerative 
cooling could satisfactorily handle was performed. The results indicated 
that 20 percent of radiation heat flux could not be handled by the available 
coolant flow but 10 percent was allowable which was approximately 654 watts/cm^ 
(4 Btu/in. sec). Therefore, the combined (convective plus radiation) maximum 
heat flux was 1153 watts/cm^ (7.05 Btu/in. ^sec) . As previously discussed, 
the carbon seeded hydrogen is capable of absorbing 95 percent of the plasma 
radiation provided certain conditions are satisfied. 

The parametric results of the 5000-kw chand>er nozzle cooling analysis are 
shown in Fig. 71. As the number of tubes (round configuration) was increased, 
the maximum tube gas-side wall temperature and the miniimim tube inside diameter 
decreased. Since a large number of tubes increases fabrication cost, a design 
with 80 tubes was selected which resulted in a satisfactory 540 K (890 F) 
maximum wall teaperature. The tube dimensions of the selected design are pre- 
sented in Fig. 72 and resulting tube wall tenperature distribution is shown 
in Fig. 73. 

In the dwsign of the Zr-Cu channel wall combustor capable of cooling the con- 
vective heat flux plus 10 percent of the radiation heat flux, the following 
guidelines were used: 
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HINIHUM TUBE INSIDE DIAMETER. INCH 




MATERIAL: 3«7 



STANCE FROM THROAT. INCHES 























1. Minimum Wall Thickness « 0.076 cm (0.03 in.) 

2. Minimum Channel Width « 0.076 cm (0.03 in.) 

3. Maximum Channel Depth- to-Width Ratio 

Hie coolant channel design for the SOOO-kw chamber is presented In Fig. 74 
^ bi-width ch«mel design (0.154 cm and 0.076 cm, or 0.06 in. and 0.03 in.) 

»as utilized to minimize coolant pressure drop. Ihls coolant circuit resulted 
in a 1.86 X 10 N/M (27 psl) coolant pressure drop and a maximum gas-side wall 
temperature of 874 K (1114 F) . The axial distribution of wall temperature, 
backwall temperature, and the coolant bulk teE,>erature is shown in Fig. 75. 

A coolant inlet pressure of 1.261 x lob N/M^ (180 psi) was necessary to pro- 
vide cooling without choked flow. A maximum coolant Mach number of 0.542 was 
obtained as illustrated in Fig. 75. 
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CONCEPTUAL DESIGN AND SCALING 


Plasma 


From the plasma core flow analysis performed by Physical Sciences Inc. (Ref. 1), 
the plasma core dimensions for a given laser intensity can be approximately 
scaled with laser power using the relationship: 

Plasma Cross-sectional Area, A a Power, P 

or that 

Plasma Diameter, D a y/l 

This relationship was used in scaling the plasma core dimensions of the 10-kw 
configuration to obtain those for the 5000-kw configuration. 

The radiation loss from the plasma is an important parameter that is not 
scalable. For the 10-kw thruster, the loss is primarily at the ultra- 
violet wavelength. However, the plasma radiation loss for the 5000-kw thruster 
would be at the ultra-violet and infrared wavelengths. Therefore, although 
the plasma formation mechanisms may be scaled with input laser power, the 
plasma radiation loss can not be scaled. 

Flow Analysis 

The relationships developed for the propellant flowrate and chamber throat 
radius (not including the boundary layer influence) result in the following 
scaling relationships with laser power level, 

<i>total “ - constant specific impulse) 

and R.J, a JT 

However, the chamber pressure, the amount of carbon seed, and the combustion 
gas temperature will also influence these parameters and must be considered 
in determining the flowrate and throat radius. 


Performance 


A comparison of the performance losses of the 10-kw and 5000-kw configura- 
tions (Table 3 and 4) reveal some scaling trends. The major influence of 
the large size difference between these two chambers is the decrease in 
boundary layer efficiency for the 10-kw thruster. The reaction kinetic 
efficiency would decrease with decrease in chamber size (thrust) for a 
fixed area ratio. However, the kinetic efficiency for the 10-kw thruster 
is higher than that of the 5000-kw thruster due to the lower area ratio. 

Optical Train 

Company- funded studies of beam propagation relating to the laser-heated 
rocket indicated that an initially annular beam will not remain annular in 
the "far field." Therefore, for space applications such as the 5000-kw con- 
figuration, the optical train configuration of the 10-kw configuration (an 
on-axis Cassegrain used as a pointer) cannot be directly used since the beam 
is no longer annular. As described in the Optical Train Section, the con- 
cepts evaluated the solid beam as well as means of forming an annular beam 
after it was received. 
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OPTICAL TRAIN ANALYSIS 


The l&ser~heeted I’ocket thruster ss configured requires s focused laser bean. 
The optical train is to receive the inconing laser bean and propagate the 
bean through a series of nirrors and/or lenses with nininun power loss and bean 
distortion to the thruster at the required focus angle. The laser bean enters 
the thruster through a solid, transparent window and focused to a spot size 
sufficient to generate the intensity required to sustain a plasna within the 
hydrogen propellant gas. 

The basic analysis for the optical train configurations (i.e., specifying 
required components; nirror design, type, cooling, mount, window/lens con- 
figuration) was made available from a con^ any -funded optical systems study 
task. 

10-kw Experimental Optical Train 

10-kw Laser Beam Characteristics . The following CO 2 laser mirror arrangement 
was assiuned to determine the laser output beam characteristics. 

Cavity length: 209 cm (82.28 in.) 

Mirror radius of curvature: 

Concave mirror: 12 m (39.372 ft.) 

Convex mirror: 8 m (26 248 ft.) 

Scraper distance from convex mirror: 13 cm (5.118 in.) 

Laser beam diameter: 5 cm (1.969 in.) 

Hole diameter: 3.3 cm (1.299 in.) 

Magnification: 1.S23 

Equivalent Fresnel Number: 3.18 

Tube Fresnel Number: 28.2 

The above configuration (Fig. 77 ) essentially represents one of the NASA-LeRC 
configurations.. Also the laser beam output was assumed to be 7,62 meters 
(25 feet) from the thruster. 



Figure 77 . 10 KW CO, Laser Configuration 






Tht* laser beam characteristics were determined using the wave XX computer 
program. This program computes modes of an unstable confocal resonator by a 
wave launching technique. The integrals required are determined using an 
average kernal technique. Using the above unstable resonator mirror con- 
figuration, a 40 X 40 grid analysis of the "far- field" was performed and then 
a finer 80 X 80 grid analysis was attempted. The analysis indicated a mode 
control problem at this low magnification of 1.523. Available computer codes 
could not propagate a laser beam with multiple modes. A minimum magnification 
of 2.0 is desirable to obtain a single transverse mode. High beam quality is 
difficult, if not impossible, to achieve with a multi-mode output beam. Two 
approaches to analyze the optical train were taken: 

1. Design an optical train for the 1.523 magnification. 

2. Design an optical train for a magnification of 2 and 
define a set of mirrors for the unstable resonator 
to give this ma^jnification. 

This latter optical train design would accommodate both magnification con- 
figurations. Since the beam output characteristics of the 1.523 magnification 
configuration cannot be defined analytically, the magnification 2 output beam 
was analyzed to establish the beam quality. It was recommended that the laser 
output beam have a magnification of 2 to ensure mode control and minimum blur 
ciicle size, and although the 1.523 magnification configuration can be used 
in the same optical train the beam quality, using the lower magnification 
particularly at focus, cannot be guaranteed. In part, the Intensity may be 
lower than that required to sustain the plasma. 

The unstable resonator mirror configuration with a magnification of 2.0 shown 
in Fig. 78, was assumed. The cavity length and the scraper distance from the 
convex mirror were the same as those for the configuration having the lower 
magnification. The CO^ laser beam output intensity and phase angle distri- 
bution for the configuration having a magnification of 2.0 are shown in Fig. 79 
and 80. 


123 






OUTPUT BEAMt^ 



Figure 80. 10 KW CO 2 Laser Rel »ttve Phase Angle Distribution (Output 

Magnification ■ 2.* } 
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10-kw Optic a l Train Confi£uratlon . As shown in Fig. 81 , several optical 
train concepts were evaluated which included: 

1. Off-axis parabola 

2. On-axis cassegrain used as a telescope 

3. On-axis cassegrain used as a pointer 

4. Scraper/ spherical concept 

Due to the problem of clearance of the flat mirror and the thruster, the off- 
axis parabola concept required a 45-degree off-axis parabola which is not 
readily available especially for high power applications. The telescope 
Cassegrain and the Scraper/ spherical concept required that the thrust chamber 
be located within the optical train. Tlierefore, the optical train concept 
selected for the 10- kw chamber was the pointer cassegrain. This concept, as 
shown in Fig. 81, utilizes a flat mirror and two spherical mirrors and tends 
to dampen any output beam tilt. 

The configuration for the 1,523 magnification output beam is presented in 
Fig. 82. Since both mirrors are spheres and since the optical system is to 
be used for on-axis object points at infinity (collimated input, parallel to 
the optical axis) the only aberration that need be considered is spherical 
aberration assuming single transverge mode operation. Therefore, a number of 
runs were made with a geometrical optics computer code that predicts the 
spherical aberration blur circle size at focus of a cassegrain system that has 
Spherical mirrors. The code predicts spherical aberration blur circle sizes 
as a function of initial ray height, initial ray slope, mirror separation, 
f/* of the secondary, and the radii of curvature of the mirrors. It was 
initially hoped that a very fast secondary mirror would give a small blur 
circle. This would allow for a very compact system. It turned out, however, 
that the system had abiiur circle size that was very sensitive to the f/* of 
the secondary (i.e., an f/1 secondary yielded blur circle diameter of 2001^ 
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when an f/2 secondary yielded one of diameter of 40y) . Therefore, an f/1.5 
secondary was used (blur circle of Sl-SSy) and the system is shown in Fig. 83. 
The blur circle will be somewhat smaller than predicted, for the following 
reason: Blur circles due to spherical aberration compare the focusing prop- 
erties of a meridional ray near the edge of the system aperture to a para- 
axial ray, which is close to the system's optical axis. In this case, however, 
the input beam has a rather large hole in it and, therefore, there are no 
rays near the optical axis. Therefore, the spherical aberration will be 
reduced because there essentially are no para-axial rays. 

Using the on-axis Cassegrain concept used for the output beam with 1.523 
magnification, an optical train was also designed for the magnification of 
2.0. As shown in Fig. 83, the on-axis spherical mirrors are larger and the 
distance between them is larger than for the configuration for the lower 
magnification. However, the configuration still lends itself to efficient 
and compact packaging. A schematic of the optical train and thrust stand is 
illustrated in Fig. 84 and 85 for both magnification values. The included 
laser beam focus angle for both configurations is 15.14 degrees. 

10-KW Chamber Laser Beam Receiving Window . A comprehensive thermal model of 
the 10-kw chamber window was constructed for the use on the Differential 
Equation Analyzer Program (DEAP) . The two-dimensional thermal model was a 
lumped parameter representation of the physical window consisting of 60 nodes, 
and capable of both transient and steady-state analyses. Each lump or node 
is connected to its adjacent nodes by the appropriate conduction, convection 
or radiation type of connector which is representative of the heat transfer 
nodes between the two locations on the physical model. Heat absorption 
arising from the penetration of a laser beam through the window is repre- 
sented by heat sources at appropriate nodes. Window back face surface nodes 
are connected to a single sink node representing the flowing hydrogen coolant. 





Heated Rocket Thruster Test Apparatus (Laser Output 













Figure 85 . jQ-KW LASER - HEATED ROCKET TEST APPARATUS 
(LASER OUTPUT BEAM MAGRIFICATIOH = 1.523) 





/ i 

The output laser beam (magnification of 2.0) was propagated through the optical * I 

train and the laser beam intensity distribution on the window is presented in | 

Fig. 86. . j 

Both zinc sulfide and zinc selenide windows were considered. Zinc selenide 

! 

has a flatter transmittance versus wave length curve at 10.6 pm. Although 

the transmittance at exactly this wave length are similar, the cost of zinc | 

I i 

sulfide is significantly lower than for zinc selenide; however the experience 
with zinc selenide led to its selection for the 10- kw chamber window material. 

Physical dimensions of the window were 3.048cm (1.200- inch) diameter and a ? 

- 1 i 

thickness of 0.6 cm (0. 236- inch) . An absorption coefficient of 0.0005cm | 

(0. 00127- inch ^) was used for the zinc selenide (ZnSe). I 

i i 

Using the sixty nodal thermal model described, the temperature distribution j 

I of the zinc sulfide or zinc selenide window was determined using window | 

intensity shown in Fig. 86. Cooling of the window by the hydrogen coolant 
flow on the chamber side and air free convection was included. The resulting 
axial and radial temperature distribtuions at 60-seconds and steady-state 
are shown in Fig. 87 and 88, respectively. A maximum temperature of 159C 
(319 F) and a maximum axial gradient of 3C (5 F) occurred on-axis. 

Beam Propagation . Laser beam propagation through the thruster ZnSe window 

resulted in the prediction that the focal point in the thruster would move 

0.18cm (0.07- inch) toward the window due to the thermal variations in the 

index of refraction of the ZnSe. The beam was propagated from the thruster 

window to focus. For the aberrated wave even in the presence of thermal phase 

distortions resulting from the window, the focal point intensity distribution 

was not appreciably different from an unaberrated wave. The on-axis intensity 

8 2 

f was computed to be 2.4 X 10 K/cm . This intensity is threv orders of 

( magnitude greater than the intensity required to sustain the plasma (3.7 X 10^ 

( 2 

j W/cm ) for this configuration. 
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Geometric optic thermal blooming analysis was performed assuming no convective 
cooling. The analysis indicated that thermal blooming was not a problem. 

Optical Train Design . The assembly drawings of the optical trains which 
are recommended for use with the 10-kw chamber are presented in Fig. 89 and 90. 

Both configurations are identical in concept and differ only in size. The flat 
mixTor and the secondary mirror are water~cooled and the primary mirror is un- 
cooled. The three mirrors are essentially contained in the same housing. The 
micrometer adjustments on the three mirrors allow tilt and the primary and 
secondary mirror tilt on two axes and a horizontal translation adjustment is 
provided on the secondary mirror. The optical train for the laser output beam 
magnification of 2.0 is shown in Fig. 89 . The one for a 1.523 magnification 
is shown in Fig. 90. The detail design drawing package composing of twenty 
drawings was prepared and was sent to the NASA-LeRC Project Manager. 

5000KW Optical Train 

Beam propagation analyses indicated that an initially annular beam will not 
remain annular in the "far field." Therefore^ for space applications such as 
the SOOO-kw configuration, the optical train configuration of the 10-kw con- 
figuration (an on-axis cassegrain used as a pointer) cannot be directly used 

since a hole no longer exists in the beam. 

A conceptual design evaluation of the orientation of the incoming laser beam 
relative to the line of thrust indicated that a 90-degree orientation pro- 
vides the optimum orientation based on beam reception from a fixed point 
(Fig. 91). Four candidate concepts employing the 90-degree orientation for the 
5000-kw optical train are presented in Fig. 92 . Two of the design concepts 
employ flat, conical and spherical mirrors, whereas the other two use aspherics 
(elliptical and parabolic mirrors). In all these concepts, the flat mirror 
could be rotated and tilted to increase the envelope of the laser beam receiver 
system. This envelope is illustrated in Fig. 93 for an attitude engine and for drag 
make- up, orbital changing, or planetary injection engines. 





Figure gg 10-kw Optical Train Assembly 
(M - 2.0) 
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Figure 90. 10-KW Optical Train Assembly (h^l.523) 
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5000- KW OPTICAL TRAIN CONCEPTS 
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10-KW EXPERIMENTAL HARDWARE DESIGN 


Water- Cooled Chamber 


Initially an evaluation of the overall fabrication technique for the 10-kw 
water-cooled chamber was conducted to develop a design incorporating forced 
convection cooling with water. The coolant passage sizes defined in the 
Cooling Analysis section were to be incorporated with a common coolant inlet 
and individual coolant outlet lines. 

In designing the water-cooled chamber several fabrication techniques currently 
used in rocket thrust chambers were evaluated and these included: 

1. OFHC copper liner with an electrode- formed nickel closeout 

2. OFHC copper liner with "0" rings within a stainless steel 
shell 

3. OFHC copper liner brazed to a stainless steel structural 
shell 


In the first fabrication technique, the electrode-formed nickel closeout 
would have provided a more than adequate structural closeout; however, for 
this particular chamber this technique would have resulted in an expensive 
solution. The second fabrication technique using "0" rings to seal the water 
resulted in a low cost technique and the capability of replaceable liners; 
however, the difficulty in sealing the electrode (igniter) region resulted in 
the elimination of this technique. A water leak into this area could result 
in a failure of the electrode mechanism to arc. The third technique of 
brazing of the copper liner and the structural steel shell eliminated the 
sealing problem with only slight increase in cost. The inner liner could still 
be replaced by machining out the "old" liner and brazing in a "new" liner. 
Therefore, this fabrication technique was selected. 
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The ignitors or plasma arc electrodes must pass through the cooled chamber 
wall; and the design of the water coolant passage arrangement in the vicinity 
of the electrodes must be such that adequate cooling be provided as well as 
sufficient space for the electrode. Therefore a small diameter electrode was 
desirable. This was also desired from a fluid dynamic aspect since an 
electrode larger in diameter than the hot-gas diameter would literally stop the 
gas flow. 

Since the plasma was predicted to exist 3 cm (1.13-inch) upstream of the 
geometric throat, water cooling was provided 5.08 cm (2-inches) upstream 
of the throat. The arc-initiating electrode mechanism mounting flanges are 
located 90 degrees from the water coolant lines. 


For testing at sea level, the 1.6:1 area ratio will result in a full flowing 

5 2 

nozzle at a chamber pressure of 3.45 X 10 N/M (50 psia). 


The chamber with the gas side wall thickness (wall thickness between the hot 
gas and the water coolant) of 0.254 cm (0.1-inch) allows the capability of 
increasing the hot gas diameter by a factor of 2. This hardware flexibility 
will allow hardware growth if the plasma size is larger than predicted. The 
assembly drawing of the 10-kw water-cooled chamber is presented in Fig. 94* 


litors (Plasma arc Initiating Electrode Mechanism) 


j The electrode design uniquely incorporated two modified TIG welding tips with 

I 0.1016 cm (0.04-inch) diameter tungsten electrodes with a micrometer 

i- 

I adjustment to set the arc gap. The basic concept is a pneumatic cylinder which 

I' 6 2 

I when pressurized (gaseous nitrogen at 3.45 X 10 N/M or 500 psig) moves the 

I two electrodes into the chamber to a preset arc distance. The instant the arc 

I occurs, the nitrogen pressure is released and the spring withdraws the electrodes 

I . flush with the chamber wall. Hydrogen pressure placed in the electrode 

I mechanism will provide positive pressure against incoming hot gas from the 
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chamber. The electrodes can be easily changed by turning the locking nut - 
by hand. Also the electrode withdrawal distance can be varied for the 
different distances required for the water-chamber and the uncooled chamber. 

The ignitor assembly is illustrated in Fig. 95. 

10-kw Heat Sink Chamber . The 10-kw heat sink chamber incorporates a thick 
copper chamber capable of 8-seconds test duration and permits a larger-than- 
predicted plasm? diameter (up to four times larger) . In addition, a graphite 
insert provides the flexibility of altering the gas contour to the contour of 
the 10-kw water-cooled chamber. The heat sink chamber utilizes the electrode 
mechanism, injector, and thrust stand designed for the water-cooled chamber. 

The assembly drawing of the heat sink chamber with graphite insert is presented 
in Fig. 96. Chamber details are shown in Fig. 97 through 99. 

Details of the graphite insert are shown in Fig. 98. Thermocouple installations 
for wall temperature monitoring are illustrated in Fig. 97 and 98. Thermocouples 
installed one-inch from the gas surface were provided at three axial locations 
(at the throat, and 2.54 and 6.35 cm upstream of the throat). A typical plasma 
viewing window installation is illustrated in Fig. 99, with a lens window 
(could be a plane window) , and a gaseous purge to prevent external hot gas 
leakage. Per a verbal agreement with the NASA-LeRC Project Manager, the plasma 
viewing window was included in the design drawings, but was not fabricated. 

Injector 

The assembly drawing of the 10-kw chamber injector is shown in Fig. 100. The 
main propellant, hydrogen, is introduced in a common manifold for the window 
coolant and the main injector. The flowrate is to be controlled with rigimesh 
orificing using 0.635 cm (1/4- inch) 100 scfm rigimesh to provide ten percent 
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flow for window cooling and the remaining hydrogen flow through the main 
injector. Also the rigimesh provides the 3.45 X IoVm^ (5 psi) injector 
pressure drop required (10 percent of chamber pressure) to maintain flow 
stability. The rigimesh was replaced with an equivalent flow rated felt 
metal due to long lead times in procuring the rigimesh. Injection velocity 
of 40 cm/sec (1.3 ft/sec) for the window coolant and 152 cm/sec (5 ft. /sec) 
for the main injection were selected to provide laminar flow across the 
window and a main propellant injection velocity closely corresponding to an 
annular velocity of the injector lip annulus (see page 77). The window 
coolant is fed throuch 72 slots 0.1016 cm X 0.0508 cm (0.04-inch by 0.02-inch) 
and the main propellant through 58 holes of 0.1189 cm (0.0468-inch) diameter. 

Thrust Stand 


The assembly drawing of the 10-kw chamber thrust is presented in Fig. 101. 

The basic concept employed in designing the thrust stand is a gravity con- 
trolled balance beam. The thrust chamber is bolted to the thrust stand 
through bolt holes in tie injector flange. The balance beam has a mechanical 
stop which prevents any damage to the flextures or load cell during mounting 
of the thrust chamber. The balance beam is essentially a channel in which 
the counterbalance weights are placed. The pivot is provided by two Bendix 
flextures. The main propellant and inlet and outlet water coolant links 
are flex- lines which originate from a main panel. The procedure for cali- 
brating the thrust stand is to: 

1. Install the pressure plate on the thrust chamber exit 

2. Pressurize main propellant to rated chamber pressure 
and water to operational pressure 

Zero out the measured thrust 
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5000-KW DESIGN CONFIGURATION 

The chamber assembly of the SOOO-Vw 

and 103 the t»o drawings illustrate T 

hydrogen flow. I„ both configurations"^"'!”^ 
combustor of ,38 varill V ““ 

he machined in a h' h i. ‘=°°*nnt channels. These coolant channels would 

elect! :l ic : “‘‘- 

reduces the gas side wall "‘‘“etion created by the channel lands which 

general the H V and the wall thermal gradients. In 

e t o h ' ie" in heavier than the con- 

ventional tubular constructinn 

channels in the high heat fl * '“"a"""'”" technique utilizing 

high heat flux, small hot-gas diameter region (subsonie to 

!d" Lri'ac?""’ -- -i-- !i::er 

thrust Chamber !thT h”'’! the best combination 

The 1 Petformance, light weight, and a high cycle life 

zrrt::::"'"'" "" " ^ - - 

Single and multi-pane windows were considered and. since the multi-pane window 

the d “‘‘'■“"tages. the single panel window was selected. Since 

the window coolant eventually flows into the chamber, the coolant p!ssu! 

US a wa>.s be higher than chamber pre.ssure 3.4S .X 10^ N/M^ (50 psia). 

ere ore. the concept of using multiple panes having decreasing coolant 
pressure could not be utilized. 

required a coolant inlet pressure of 1.242 X 10^ N7M“ 

(80 psia). necessary to maintain reasonable coolant velocities (M=0 54^) 

This results in a 1,0.^6 X 10^ ’ ^'•^>4.). 


N/M‘ 


(150 psia) coolant pressure which must be 


regulated through a series of orifirt^c 

T i- Y ifiS y,.,? . iTices prior to injection into the 

* ‘ ‘ (^*^ psia), chamber pre.ssure. 
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In the first configuration (Fig.i 02 ) the coolant flow exiting the cooling 
jacket splits 60 percent to the carbon seeder and the remaining 40-percent to 
the main hydrogen injection ring (30-perc.2nt) and the window cooling (10 per- 
cent). In this configuration, the carbon .eeded flow is "sandwiched" between 
streams of hydrogen flow at the main injector. 

As discussed in the plasma analysis section, to absorb a major portion of 
the plasma radiation, the flow velocity of the carbon seeded layer must be low 
compared to the main propellant flow velocity. Therefore, a second chamber 
configuration utilizing essentially two axial locations for the injection of 
the seeded flow was designed as shown in Fig. 103 . The seeded flow injection 
ports consist of holes through the land of the coolant channels . This con- 
figuration would enable better control of the cooling layer velocity in the 
region where the plasma exists. 
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Figure 103 5000*kw Chamber Configuration 
(single Carbon Seeded Flow 
Injection) 
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TASK II: FABRICATION 


10-KW THRUSTER 

All components of the 10-kw thruster were fabricated using conventional 
manufacturing techniques currently using in fabricating rocket engines. 

Water-Cooled Chamber 

The water-coolod thruster component (copper liner, structural shell, and 
ignitor flanges) machining was performed by Tri-Models and the brazing and 
welding were performed at Rocketdyne. 

The OFHC copper liner, the stainless structural shell, and the ignitor 
flanges were machined in the first stage of the water-cooled chamber fab- 
rication. These components are shown in Fig 104. In the first braze cycle 
the copper liner and the outer structural shell were brazed together. Sub- 
sequent to this braze cycle the chamber coolant channels were leak tested. 

This test revealed only minor inter-passage leakage and no leakage into the 
electrode region and no external leakage. In the next step of the fabrication 
process, the ignitor flanges were brazed to the chamber assembly. Next the 
coolant inlet manifold cover plate and inlet line were TIG orazed to the 
chamber. The internal chamber contour was then machined and the ignitor ports 
final machined. Next the nine coolant outlet lines were welded to the chamber 

assembly. 

A thermocov 5 >le was installed approximately 6.35 cm (2.5 inch) from the chamber 
throat and 0.508 cm (0.20 inch) from the gas-side surface in the uncooled 
section of the chamber. This thermocouple will provide an indication of 
Chamber overheating in this region. 

Upon completion of the thruster assembly (Fig.i05) coolant channels, inlet 
manifold, and coolant outlet lines were successfully hydrostatically tested 
to 3.45 X 10^ N/M^ (500 psig) without leakage. 
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gure 105 lO-kw Water-cooled Chamber 




Once the ignitors were assembled the fit with chamber ignitor ports were 
checked and the ports were honed to provide a correct fit. The complete 
water-cooled chamber assembly with injector and ignitors is shown in Fig. 106. 

Ignitors (Plasma ARC Electrode Mechanism) 

The ignitor components which included the ignitor housing, electrode key, 
retracting spring, conductor rod and insulator, nut wrench, micrometer housing 
sleeve, and piston (Fig. 107) were fabricated by Acro-Tech. These components 
were assembled at Rocketdyne and the one of the final assemblies is shown 
in Fig. 108. The electrode locking mechanism was checked functionally to lock 
and retain the tungsten electrode and unlock. The pneumatic piston actuation 
of the ignitors were checked using approximate 1.034 X 10° ti/VT (150 PSIG) 

Argon. Although the specified actuation pressure was 3.45 X 10° N/M (500 PSIA) 
gaseous nitrogen, these actuation tests indicated that a lower gaseous nitrogen 
perhaps 1.379 X 10^ to 2.068 X 10^ N/M^ (200 to 300 PSIG) would be sufficient. 

A lower gaseous nitrogen pressure would reduce the electrode retracting time. 

Uncooled Chamber 


The uncooled OFHC copper chamber and the ATJS graphite insert were fabricated 
by Tri-Models. The uncooled chamber was machined from a solid copper billet. 
The ten thermocouples on the uncooled chamber were installed at Rocketdyne. 

The ignitors, spacers and chamber were assembled for an assembly check as 
shown in Fig. 109. The uncooled chamber assembly complete with Ignitors and 
injector is shown in Fig. 110. The cooled and uncooled chamber injector orifice 
plates are shown in Fig. 110. 

Injector 


The injector components (Fig. ill) which include the injector body, orifice 
plates, felt metal ring, window seals, window coolant ring, and window cover 
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plate were fabricated by Acro-Tech. Assembly of the injector and installation 
of the chamber pressure measuring tap was performed at Rocketdyne. The in- 
jector assembly is shown in Fig. 110. 


10- KW THRUST STAND 

The 10-kw thrust stand components including the balance beam, base, interface 
panel, pitot clevis, and thruster pressure plate were fabricated by Tri-Models. 
A check on the chamber fit on the thrust stand indicated that the thrust 
stand was not to print and therefore required a modification so that the 
chamber would fit. The thrust stand assembly is shown in Fig. 112. 
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